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Abstract 

In this article, we develop a general CARA framework to study optimal execution 
and to price block trades. We prove existence and regularity results for optimal 
liquidation strategies and we provide several differential characterizations. We also 
give two different proofs that the usual restriction to deterministic liquidation strate- 
gies is optimal. In addition, we focus on the important topic of block trade pricing 
and we therefore give a price to financial (il) liquidity. In particular, we provide a 
closed-form formula for the price a block trade when there is no time constraint to 
liquidate, and a differential characterization in the time-constrained case. Numerical 
methods are eventually discussed. 

Introduction 

A general issue for stock traders consists in buying or selling large quantities of 
shares within a certain time. Unlike for small trades, a trader executing large blocks 
of shares cannot ignore the significant impact his orders have on the market. A 
trader willing to sell instantaneously a large quantity of shares would indeed incur 
a very high execution cost or could even be prevented to succeed in selling because 
of limited available liquidity. 
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In a nutshell, traders face a trade-off between price risk on the one hand and both 
execution cost and market impact on the other hand. Traders liquidating too fast 
incur high execution costs but being too slow exposes to possible adverse price fluc- 
tuations, effectively leading to liquidation at lower-than-expected prices. For that 
reason, traders usually split their large orders into smaller ones to be executed pro- 
gressively within a certain time window. Research on optimal execution^ has initially 
been dedicated to this issue of optimally splitting those large orders|§ 

The above trade-off between execution cost, market impact and price risk first ap- 
peared in the economic and financial literature with Grinold and Kahn [TS] and has 
been widely studied since then thanks to the framework developed by Almgren and 
Chriss in their two seminal papers [21 E]. It is noteworthy that the literature on op- 
timal liquidation started beforehand with Bertsimas and Lo [TJ, but their research 
focused on the minimization of expected execution costs only, consequently ignoring 
price risk that plays an important part in practice. 

In the last ten years, the framework proposed by Almgren and Chriss has largely 
been used in practice and generalized either to better fit real market conditions or 
to enlarge the scope of modeling possibilities. Initially developed in discrete time 
with linear execution costs and within a Bachelier model for the price, it has also 
been considered in continuous time and generalized to allow for nonlinear execution 
costs and random execution costs [5] J§ Black-Scholes dynamics for the price has also 
been considered and attempts to generalize the model in other directions have also 
been made, for instance to take account of stochastic volatility and liquidity [TJ. 
Discussions on the optimization criterions and their consequences on optimal strate- 
gies are also very present in the literature (see for instance pE], [H], [20] and [26]). 
In this paper, we consider a Von Neumann-Morgenstern expected utility framework 
and the specific case of an investor with constant absolute risk aversion^ This 
CARA framework has been studied in [21] in which the author proved that going 



1 Optimal execution may sometimes be referred to as best execution or optimal liquidation. In 
this paper, we only focus on liquidation but the issues of buying large quantities of shares can be 
tackled in the same way. 

2 Today, new strands of academic research have developed. Following the seminal paper by 
Obizhaeva and Wang |21) , many authors model the dynamics of the order book instead of having 
a statistical view on execution costs. Also, the focus of research has slightly moved from time 
scheduling to the actual way to proceed with execution. Liquidation with limit orders - see [H I16[ 
[17] - and dark pools - see [HI [19] - are now important topics. 

3 The proofs in [5] are however limited to the case of a C 2 execution cost function and C 2 
trajectories... two unrealistic assumptions as we shall see in the text. 

For short periods of time there is no real difference between Bachelier dynamics and Black- 
Scholes dynamics. 

5 Very interesting results in the case of IARA and DARA utility functions are presented in [23] , 
For practical applications, choosing the appropriate level of risk aversion is already a complex task 
and we shall not go beyond CARA utility functions. 

6 We repeat below their proof that deterministic strategies are optimal. Another proof is also 
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from deterministic liquidation strategies to stochastic liquidation strategies does not 
bring any advantage. 

In this paper, we develop a framework that is general enough to cover most nonlin- 
ear execution cost functions used in practice. In addition to nonlinearity, we cover 
the case of (deterministic) time-dependent execution costs ; an important case to 
account for the shape of market volume curvesJll However, as far as price dynamics 
is concerned, we focused on the case of a Bachelier model with a drift given by a lin- 
ear permanent market impactfl Optimal liquidation strategies are discussed, along 
with the important problem of block trade pricing. The academic literature indeed 
focuses on the rhythm at which liquidation should be carried out but ignores the 
global cost of liquidation. We regard this issue as an important one, both to know at 
what price to buy a large portfolio and for risk managers to evaluate liquidity premia. 

We present the general framework in Section 1. Section 2 presents existence, unique- 
ness and regularity results for the optimal deterministic liquidation strategy. Dif- 
ferential characterizations are provided that are useful when it comes to numerics. 
In addition, we provide a general proof of the reduction of the problem to a first 
order ODE when market volume is constant. Section 3 gives a first proof that no 
stochastic liquidation strategy can improve the best deterministic strategy. Section 
4 is dedicated to block trade pricing. We first prove regularity results for the value 
function^ We then determine the price of a block of go shares given market condi- 
tions and the fact that the shares have to be liquidated within a certain time window. 
The partial differential equation associated to the pricing problem (and to the initial 
liquidation problem) has a singularity at final time and we provide a way to approx- 
imate the price of a block trade using a relaxed problem that does not have any 
singularity. Interestingly, we also provide a closed-form expression for the price of a 
block trade when no time window is imposed for liquidation afterwards. Section 5 
uses the relaxed problem of Section 4 in which liquidation is not enforced to provide 
another proof that no stochastic strategy can do better than the best deterministic 
strategy. Finally, Section 6 is dedicated to a discussion about numerical methods to 
find the optimal liquidation strategy and/or the price of a block trade. Numerical 
methods are rarely discussed in the literature (see [12] for one of the rare examples) 
but we shall see that some natural discretizations may not provide satisfying results. 

provided in Section 5 using a completely different approach. 

7 The available liquidity is not the same at each hour of the day. In Europe, in addition to the 
importance of the near-opening and near-closing periods, the opening of the New York marketplace 
is an important event. 

8 The choice of a Bachelier model over a Black-Scholes model is mainly dictated by mathematical 
considerations. However, there is no real difference from the modeling side because liquidation 
problems are short-term problems. 

9 In [53], the authors discuss the properties of the value function. However, contrary to what 
they claim, the only thing they prove is almost everywhere differentiability and that the Hamilton- 
Jacobi is solved almost everywhere. 
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1 Setup and notations 

A trader with a portfolio containing^ go > shares of a given stock is willing to 
liquidate his portfolio over a time window [0, T]. 

Let us fix a probability space (ft, J 7 , P) equipped with a filtration (J : t)telo > T\ satis- 
fying the usual conditions. We assume that all stochastic processes are defined on 

(7- t ) t6[0l T],P). 

We introduce the set V(t,T) of progressively measurable processes defined on [t,T] 
and the set A = {(u t ) te[0 ,T] G V{0,T), (t,u) \-> f*v s (u)ds G L°°([0,T] x For a 
control process (ft)te[o,T] £ -A - representing the velocity at which the trader sells 
his shares -, we denote by (qt)te[o,T] the number of shares in the portfolio where q t 
is given by: 

Qt = qo~ / v s ds. 
Jo 

The trades impact the market price in two distinct ways. Firstly, there is a per- 
manent market impact (assumed to be lineai0) that imposes a drift to the price 
process (S t )te[o,T]- 

dS t = adW t - kv t dt, a > 0, k > 0. 

Secondly, the price obtained by the trader at time t is not St because of execution 
costs. To model these execution costs we introduce a function L G C 1 (R, R) verifying 
the following hypotheses!^! 



L(0) = 0, 

L is an even function, 
L is increasing on M + , 
L is strictly convex, 

3e > 0,3A,a > 0,B,b> 0,Vp G R,L{p) > A\p\ l+e - B and \L'(p)\ < a\p\ e + b. 



10 The case qo < can be treated using the same tools. 
"See [TJ. 

12t 



2 We want to cover the cases L(p) = r]\p\ 1+ ^ for rj > and <fi > 0. In particular, the function 
L is not C 2 when (f> < 1 (a common case when the function L is calibrated on real data), and 
although it is a C 2 function when <fi > 1, it does not verify in this latter case that L" > 0. 
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We also introduce a market volume process (Vt)te[o,T\ assumed to be continuous, 
deterministic and such that 3V > 0,V > 0, Vt G [0,T\,V< V t < V. 

This allows, for v G A, to define the cash process^] (Xt)te[o,T\ as: 



where the execution cost is divided into two parts: a linear part which represents 
a fixed cost (if; > 0) per share - linked to the bid-ask spread and the tick size -, 
and a strictly convex and superlinear part modeled by L, to account for the cost of 
trading fast. 

We define our objective function for v G A as 



where 7 > is the absolute risk aversion parameter of the trader. 

The goal of this paper is to solve two problems: 

• We want to find an optimal strategy to liquidate the portfolio within a certain 



time T. In other words, we want to maximize J over Aq = < v G A, J Q v s ds = q 



• We want to give a price to a block trade. In other words we want to determine 
the maximal price an agent should buy a portfolio of go shares of the stock, 
given that this portfolio will have to be liquidated on the market. 

2 Optimal trading curves 
2.1 Introduction 

In this section, we are going to consider the case of deterministic strategies for op- 
timal liquidation. Deterministic liquidation strategies are simpler to analyze than 
general stochastic ones because the cash process at the terminal time T (that is 
Xt) is normally distributed when the liquidation strategy (v t )te[o,T] is deterministic. 
As we shall see in Section 3, it turns out that no stochastic strategy can improve 
the best deterministic strategy and that the restriction to deterministic strategies is 
therefore not a real restriction. 

To start, let us define the set *4det of deterministic strategies in A and the set v4o,det 
of deterministic (liquidation) strategies in Aq. 



The cash process may be equal to —00. 




J(v) = E[-exp(-7Xr)] 
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Depending on whether v is in A, Adet or *4o,det, we have different expression and 
properties for the cash process. The following proposition and its corollary state 
these properties. 

Proposition 1. Let us consider v G A andt G [0, T\. We have: 

X t + 1tS t - ^q 2 = QoSo - 7^0 ~ j V sL ^y) ds-ip j \v s \ds + J aq s dW s . 
In particular, if ' v G Adet, if Jq V s L (y^j ds < +oo, and if J* \v s \ds < +oo then, 

X t + q t S t - ^q 2 ~ M \q S - |<?o ~ J V s L ds - ip j \v s \ds, a 2 J q 2 s ds^j . 

Proof: 

By definition: 

X t = j v s S s ds - j V S L \y)ds-i) / \v s \ds 



JO \ y s/ JO 



t rt rt / X rt 



qoS -q t S t -k / v s q s ds+ I aq s dW s - J V S L ^— J ds - ip J \v s \ds 

"t /„. \ rt rt 



= -q t S t + |g 4 2 + go^o - - J V S L ds - ip \v s \ds + aq s dW s . 

When (g s )se[o,t] is deterministic, if JqV s L (jM ds < +oo and J* \v s \ds < +oo, then 
X t + qtSt — fgf is normally distributed with mean go^o — \q 2 — JqV s L (j^J ds — 
tp J* \v s \ds and variance a 2 J* q 2 ds. □ 

A straightforward corollary of Proposition [I] is that Xt is normally distributed in 
the deterministic case when q? = 0. More precisely: 

Corollary 1 (Distribution of Xt)- Let us consider v G Ao,det- If JqVsL yyj ds < 
+oo, and J* \v s \ds < +oo, then: 



X T ~ Af (q S - k -ql - £ V S L f^r) ds 



T 

ip I \v 8 \ds, a 2 / q 2 ds ) . 



Using the above corollary and the Laplace transform of a normal distribution, we 
can find a closed-form expression for the objective function J: 
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Corollary 2. Let us consider v G Ao t det- Then: 

J( v ) = - exp ^-7 (q S - 7^0 ~ V s L (yj ds ~ ^ f Q \ Vg \ ds ~ ^ 7<j2 J Q q * ds 

2.2 Existence, uniqueness and characterization of a mini- 
mizer 

We can then define a new objective function for v G »4o,det by 

I (V)= J V * L (yj d s + i> J \vs\ds + )pa 2 J q 2 s ds, 

so that: 

J( v ) = - exp ^-7 (q S - 7^0 ~ H v ] 

Then 

sup J(y) = - exp ( -7 ( g So - 77<?o ~ inf 
and the maximizers of J correspond to the minimizers of I. 

We are going to study the function I but instead of regarding I as a function of v, 
we regard I as a function of g. More precisely, we define: 







(V S L (^)+^|g( S )| + h«-V( S ))^, 



where the Sobolev space W^+^i, £2) is the set of functions q G W l,1+t (ti,t 2 ) sat- 
isfying g(ii) = ai and qfa) = a>i- It is a natural space for our problem because of 
the assumptions on L. 



The following theorem states that there exists a unique minimizer of 1^ in W ' (0> T) 

Theorem 1 (Existence and uniqueness of a minimizer) . There exists a unique min- 
imizer q* G W^ o +€ (0,T) of the function 2^. 77ms minimizer is a nonnegative and 
decreasing function. 

Proof: 

Let us first consider the function qi ine G W^ +e (0, T) defined by qu ne (t) = q (l — |r) . 
2^/) (ftme) < +°° an d therefore 'vai w i,i+ e ^ QT ^ 2^ < +oo. 

Now, let us consider a sequence (q n )nen of functions in W^'^^O, T) such that 
lini n -*-+oo 2^(g n ) = inf lv -i,i+ef T}-^/' - 
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For n G N, we can define q n by q n (t) = inf s < t (g n (s)) + . We have q n G W q ^ e (0,T) 
with q n = q n ^{q n >o,q n =q n }- Hence, using the monotonicity property of L, we obtain 
that T$(q n ) < 2^,(<Zn) an d we therefore have that lim n _ 5>+00 I i /,(q ; ri ) = inf^i.i+e^^X^,. 



Now, we have: 

fT 



A 



- BV, ds 



Hence sup n6N H^nllLi+^o.r) < +00. Now, since \q n \ < q , we know that (q n ) nt =n is a 
bounded sequence of Wg^~ e (0,T). Hence, there exists q G Wg^ o +e (0,T) such that, 
up to a subsequence, q n — > q in C°([0, T]) and q n — ^ g in L 1+e (0,T). 



Now, 



y(M) 



< a 



+ 6 so that the function s H> V [ ) is in L 1+ t(0,T). 



Hence, ^ L' (<?n( s ) —q{s))ds — >■ and we have using the convexity of L that: 



2y>(<?) 



< 



' / V x L(^)+i>\q( S )\ + l 1 a 2 q 2 ( S ))ds 



T 



V. 



\ VL ( ] c/.s 



1 2 



q 2 (s)ds 



[Q n (s) ~ q(s))ds 



1 /" T 

+ ^% + ip° 2 J q 2 (s)ds. 



Hence: 

Z^g) < liminf f V S L ds + i>q Q + Ua 2 f q 2 {s)ds 

< liminf/ f %^^) cis + ^go + ^ 2 lim / q 2 n (s)ds 

< liminf Z^(q n ) 

< inf . 

< !o +E (0>T) 



Hence g* = q is a minimizer of 2^ and, as announced, it is a nonnegative and de- 
creasing function (as the limit of a sequence of such functions). 

Coming to uniqueness, if there were two different minimizers, q* and q**, then we 
would have 2^ ( 9 + £ ) < \X^{q*) + |2^(g**) = Z$(q*) because L is convex and 
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q i— > q 2 is strictly convex, in contradiction with the optimality of q*. Hence, q* is 
unique. □ 



Now, let us notice that if q G W^'^^O, T) is decreasing, then: 



'0 

"T 



'0 



• ^ cruris) . 



V. / 2 



2b (?) + V>?o- 



Hence, g* is the unique minimizer of X and it does not depend on ip. From now, we 
shall refer to X as X and q* will be the unique minimizer of X. 

Using this property of q*, we are going to derive two different differential character- 
izations of q* and eventually prove that q* G C 1 ([0,T]). 

We first exhibit the Euler-Lagrange equation for q*: 

Proposition 2 (Euler-Lagrange equation for q*). q* is the unique weak solution of 
the differential equation: 

iL'( q ^-)=^q(t), q(0) = q , q(T) = 0. 



dt \V t j 
Proof: 

For <f) G C^fO, T]) with 0(0) = 0(T) = 0, we introduce :eGl4 X(g + e0). 
The optimality of g* implies that i' ? *^(0) = 0, that is: 

/ T ( L ' (If) ' (t) + 7(T2?(t)0(t) ) rft = °" 
This is exactly the weak formulation of the Euler-Lagrange equation. 

Now, if q is a weak solution of the Euler-Lagrange equation, because i is convex, we 
have: 

Hence, q = q* ■ □ 
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Associated to this Euler-Lagrange characterization, we exhibit the Hamiltonian char- 
acterization of q*. For that purpose, we introduce H, the Legendre transform of L: 
H{p) = sup pgR pp — L(p), and we recall that the strict convexity of L implies that 
H is a C 1 function. 

Proposition 3 (Hamiltonian equations for q*). Let us introduce p*(t) = V ■ 
Then, (p*,q*) solves the following system of equations: 

p{t) = -fo- 2 q(t) 

m = vyrw)) 9(0) = *- q< ) = °' 

Moreover, if there exists a couple (p, q) of two C l functions satisfying the above 
equations, then q = q* . 

Proof: 

By duality, p*(t) = L> (^f) =► = H'(p*(t)) =► q*(t) = V t H'(p*(t)). 

Now, if (p, q) is a couple of C l functions satisfying the above equations, then, by 
duality, p(t) = V (®f). Hence: 

i L '( q 4r)=^Mt), q(0) = q o , q(T)=0, 



dt \ V t j 

and the conclusion follows by the above Proposition. □ 

A consequence of this differential characterization is the following: 
Corollary 3. q* G C\[0,T]). 

We shall see in the next subsection that there exists cases where q* ^ C 2 ([0, T]). This 
means that one has to be careful when considering the Euler-Lagrange equation. 
Now, since q* is C 1 , v* = —q* G A),dct, we can gather all the results to obtain a 
theorem about the minimizers of J: 

Theorem 2 (Existence and uniqueness of a minimizer for J in «4o,det)- There exists 
a unique minimizer v* of J in Ao,det ■ 

Moreover, t (->■ v% is a continuous and nonnegative function. 

2.3 The special case of a flat market volume curve 

In this subsection, we consider the special case V t = V. In that case, and under 
smoothness assumptions on L, the Euler-Lagrange equation can be simplified and 
boils down to a first order ODE. 
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Theorem 3 (First order ODE for the optimal liquidation strategy) . Let us assume 
that V t = V. Moreover, let us supposS that L G C 2 (R*) with Vp G R*, L"(p) > 0. 
Let us introduce the strictly increasing function r : p G R+ h-> pL'(p) — L(p). 

m If £ i— )■ pi7|2) ^ s no ^ integrable in £ = i/ien t/iere exists C > s«c/i g* zs 
uniquely characterized by: 



q*(t) = -Vr~ l 



7<7 
2V 7 



g*(t) 2 + C 



with <f(T) = 0. 



//£ H- r _ 1 1 ^ 2 ^ ) is integrable in £ = 0, t/ien ie£ us define gQ im 



g im = +00, 



Jo 



_lim . 

yo 1 



1 / T CT c2 



d£<T 



-<i£ = T, otherwise. 



For q > q l Q m , there exists C > suc/j i/iat g* is uniquely characterized 
by: 



q*(t) = -Vr 



-1 



2V 



q*(t) 2 + C 



with g*(T) = 0. 



For qo < g im , g* is uniquely characterized by: 



q *(t) = -Vr" 1 



2^ 



with g*(0) = g . 



Moreover, if q G C 1 ([0,T]) wif/i g(0) = g ariC ^ 9 CO = solves, for some constant 
C > 0, the equation: 

q(t) = -Vr- 1 



10 q (t) 2 + c 



2V 



then q = q* . 
Proof: 



The function r : p G R+ (->■ pL'(p) — L(p) verifies: 

• r(0) = 0, 

• r is strictly increasing (because L is strictly convex), 

• r is continuous, 



14 This hypothesis is verified by the execution cost functions we mentioned in the introduction. 
In particular, it is important to exclude p = in the assumption because L(p) = r]\p\ 1+ ^ is not C 2 
in 0. 
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• r is C l on (0, +00), with r'(p) = pL"(p), 

• lim p ^ +00 r(p) = +oo@ 

• r is a bijection from K + to R+, 

• r^ 1 is continuous, 

• r _1 is C 1 on (0,+oo). 

Now, let us consider a function g G C 1 ([0, T]) with g(0) = go, q{T) = and satisfying 
for some constant C > 0: 



7<r z 
2V 



q(tf + G 



We argue that q = q* . To prove this let us consider two cases: 
Case (i): Vt < T,q(t) > 0. 

In that case, q is C 2 on [0, T), g < on [0, T), and we can differentiate the equation 
r = ^q{t) 2 + C, to obtain for t G [0,T): 



V \ V 

q(t)q(t) L „{ <i(t) 



70" 



q(t)q(t) 



V 2 



V 



V 



q(t)q(t). 



Dividing by 4^, we get: 



i.e.: 



V 



dt 



V 

m 

V 



7 a 2 g(t), 



icr 2 q(t). 



Because L is even, this is equivalent to the Euler-Lagrange equation: 

d r , ( q(t)' 



u 

dt \ V 



ya 2 q(t). 



We then conclude that q = q*. 



Case (ii): 3r < T, g(r) = and Vt < r, g(t) > 0. 



15 If indeed r was bounded then Vp > 1, (j^jf^J 
L(l) + 1 1 r I |oo, contrary to our assumptions on L. 



_ r (p) 



— 



Hence, Vp > 1, ^ < 



12 



In that case the first thing to notice is that necessarily C = and q(t) = 0, V£ G [r, T]. 
We know indeed that q is decreasing (with q < — Vr~ x (C)) and that q(T) = so 
that g must remain equal to zero once it touches 0. 



Now, we can apply the reasoning we applied for case (i) to the interval [0, r). Hence, 
V* G [0,r): 

Now, because q = on (r, T] , the above equation also holds on (r, T] . Consequently, 
q is a weak solution to the Euler-Lagrange equation and therefore q = q* ■ 

Now, we need to prove that Vg > 0, 3C > 0, 3q G C^QO,^) with q(0) = q and 
q(T)=0 such that: 



q(t) = -Vr- 1 



IV 



q(tf + C 



For that purpose, let us introduce the following backward Cauchy problem in [0,T] 
for a given constant C: 



(S c ) q(t) = -Vr' 1 



2V 



q (ty + C 



with q(T) = 



t i — y q(t) = is a solution of the equation (£ ). However, this solution is not unique 
whenever £ h-> is integrable in £ = 0. 

If ^ i — >- r _ 1 1 ^ 2 ^ is integrable in £ = 0, let us implicitly define the function g hm by: 



I 



dq 



Vr-* (^g 2 ) 



T-t, 



forte T — 



Vr~ 



T 



This function defines a solution to (So). 



Jo 



Vr' 



\ < T, then for all positive q , we can consider the equation 
K) m = -Vr- 1 {^q(t)^ with q(0) = q 



We know that £ > t— > —Vr J * s continuous and decreasing. Therefore, 

using Cauchy-Peano theorem for existence and monotonicity for uniqueness, there 
exists a unique solution q to (£' qo ) on some interval [0,£ 90 [. Because this solution is 
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a decreasing function, it either merges with q hm at some pointl 16 ! t < T and is equal 
to at time T, or it reaches for some t G [0,T] (and stays equal to after time 
t). Hence, the solution q of (£ ' ) is a solution to (£q) on [0, T] and solves our problem. 

POO 

Now, if / p-^x — v- > T then q hm is defined on [0, T] and the above reasoning 

Jo v^i^i 2 ) 
applies but only for q Q < g lim (0). 

The problem for a given q > g hm (0) is then to find C > so that a solution of (8c) 
reaches qo in t = 0. To this purpose, we define, for C > 0, the following family of 
equations: 

VoJ m = -Vr- l (^q{t) 2 + C^ with q(0) = q 

Let us consider for a given C > the maximal solution of (8 Cqo ) (the Cauchy 

problem is well-posed because f i-> -Vr- 1 (f^£ 2 + is locally Lipschitz). Since 

g(t) < — W -1 (C) < 0, we know that there exists tc so that q(tc) = 0. Now, since 
q is decreasing, we can write 

f qo dq 
tc ' 



rgo 
Jo 



'o Vr- 1 (f^g 2 + C 
C i — y tc is continuous with: 



lim t c = / / : x > / / : x = T, 



C->0 J y r -i{^ q 2\ J Q Vr -l(lg q 2 



and 



lim tc = 0. 



Hence, there exists C > so that = 7\ For this C, the solution q of (8 c>qo ) solves 
(£ c ) with q(0) = q . 

Our problem is then solved in the case where £ i— >■ pr?|2) is integrable in £ = 0. 

Now, coming to the case where £ H- pijb^ is not integrable in £ = 0, we can proceed 
as above. If we indeed fix q > and consider the family of equations ), we 
know that for any C > there exists tc so that the solution q of (8 C ) verifies 

g(tc) = o. 



16 In fact it cannot merge at time t < T because the backward Cauchy problem associated to 
the junction point would be well-posed since £ H> —Vr^ 1 (^r^ 2 ^j ^ s locally Lipschitz on M!j_. 
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Since q is decreasing we can write 

" qo dq 



t 



c 



C t-7- t c is continuous. Since £ i— > r -i^ is not integrable in £ = we have 
linic^o tc = +oo. Also, as above, lim<7_s. +00 tc = 0. 

Hence, there exists C > so that tc = T. For this C, the solution q of (£p ) solves 

(£c) with g(0) = g - □ 

This simple first order ODE can be used to find in closed-form the optimal liquidation 
strategy in at least two cases. The first case corresponds to the cost function L(p) = 
rjp 2 introduced by Almgren and Chriss [3j. In that case we have the following 
proposition: 

Proposition 4 (The Almgren-Chriss case). Let us consider L(p) = rjp 2 . IfVt — V 
then q* is given by: 

smh(^(T- t)) 

Proof: 



2 



This is a pure application of the preceding theorem with C = l&r 1 f y - D 



su,l,-'| J^fT 



The second case corresponds to L(p) = rjp < , (5 > 0), for q small. 

Proposition 5 (Super-quadratic execution costs for small inventories). Let us con- 
sider L(p) = r]p 2+s for 5 > 0. Let us assume that V t = V . 

Then, for q < (4*)^ T^V^ (^5) ) *> <?* 18 9 men h V : 



2+5 
S 



Proof: 



This is a straightforward application of the preceding theorem with C = 0, the con- 
stant g im being equal to (^)^ T^V^ (^))"- □ 

The latter case exemplifies the limited smoothness of q* since q* is C 1 but not C 2 . 
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3 Deterministic strategies versus stochastic strate- 
gies 



In the above section, we only considered deterministic strategies. We now provide 
a proof that no stochastic strategy can improve the best deterministic one. This 
simple proof is based on the use of Girsanov Theorem and has been proposed in 
|24J. We repeat the argument here for the sake of completeness. 

Theorem 4 (Optimality of deterministic strategies). 

sup E [— exp (—7X7)] = sup E [— exp (—7X7-)] . 

Proof: 



For any v G Aq, we know that 



X T = q S - ^ql - J V S L f^j ds - ip J \v s \ds + aq s dW s . 



Hence: 



E [- exp (-7X T )] = - exp ( -7 ( q S - -q 



xE 



exp l 7 



V,L 



Vs 



ds + ip \v s \ds I I exp ( — 70" / q s dW : 



Hence, if we introduce the probability measure Q defined by the Radon-Nikodym 
derivative^] 

1 



dF 



exp I -70- / q s dW s - --y 2 a 2 I q A s ds 



then: 



E [- exp (-7X T )] = - exp ( -7 ( q S - -q 



xE^ 
Hence: 



exp 7 



ds + i/j I \v s \ds ) ) exp ( 7p 2 °~ 2 \ it^s 



E[-exp(- 7 X T )] = E c 



Now, let us fix u G Q. If t i-> q t (u) G W 1A+e , then X^{q{u>)) > X^(q*). Otherwise, 
Ijjj(q(u))) = +00 and we have X^(q(u)) > X^(q*). 
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We can apply Girsanov theorem since q is bounded. 
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This leads to 



E [- exp (-7-Xr)] < - exp ( -7 ^q S - t^o ~ ^(Q*) 

i.e.: 

E [- exp (-jX T )] < sup E [- exp (-~/X T )] . 

We then obtain 

sup E [— exp (—■yX T )] < sup E [— exp (—jX T )] . 

Since the converse inequality holds, the result is proved. □ 

This theorem states that no stochastic liquidation strategy can do better than v*, 
not only in Ao^t but more generally in Aq. Hence, we solved our first problem 
and we will present, in Section 6, numerical methods to approximate the optimal 
liquidation strategy. 

We will also present another proof of the result of Theorem H] based on a completely 
different idea in Section 5, more precisely using a relaxed problem to remove the 
singularity 




4 Block trade pricing 

In addition to the optimal liquidation strategy, we are interested in the determination 
of a price for a block trade. Given the state of the market, we want to determine 
the maximum price at which an agent would be ready to buy a portfolio containing 
go shares when he is forced to liquidate the portfolio on the market within a certain 
time T. This maximum price P(T, go , So), in our expected utility framework, is 
obtained by the indifference pricing approach: 

P(T,q ,S ) = --log ( - sup J(v) J = q S - -ql - ipq - 1{q*) , 

where q* is the optimal trajectory associated to the liquidation of a portfolio with 
go shares on the interval [0, T], as in Section 2. 

In practice, one can first compute the optimal liquidation strategy and then com- 
pute the price P(T,q , S ). This is a two-step approach. Another approach consists 
in computing directly the price of a block trade. This one-step approach (or direct 
approach) is based on the value function of the optimal control problem of Section 
2 and relies on the theory of first order Hamilton- Jacobi equations. 
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We introduce the value function of the control problem of Section 2: 



Or(t,q)= inf / (V 8 L ^ + - 7 aV(«) ) da. 

Since we bounded ourself earlier to the case of positive inventories, this function is 
a priori defined on [0, T) x . We rather consider it on [0, T) x R - and we notice 
that Vt G [0, T), fl T (t, 0) = and V(t, g) G [0, T) x R, 6 T (i, -q) = 6 T (i, q), because L 
is even. 

The introduction of 9t allows to write P(T, go, •S'o) as go'S'o — fgo — ^Qo — #t(0, go). 
More generally, we can write the price P(t, q, S) to liquidate q shares on a time 
window of length t when the price of the stock is S as 

P{t, q, S) = qS- \q 2 -^q- 9 T (T -t,q),VT>t 

We then see that the study of block trade pricing boils down to the study of the 
value functions fly's. 

4.1 Study of 6 T 

We are now going to study the properties of fly. For that purpose, let us denote 
q^ G W-j' +e (t, T) the function at which the infimum in the definition of fly(t, q) is 
attained. Such a function does exist by Theorem [T] since we can replace the couple 
(0, g ) by the couple (t, q) - and it is a C 1 function on [t, T\. 

A first important result is the behavior of fly(t, g) when t — v T. Proposition Instates 
indeed that fly has a singularity at time T, because of the liquidation constraint 
q T = 0. 

Proposition 6 (Singularity of fly). 

t-*T [+oo, otherwise. 

Proof: 

fly(t, 0) = 0. The result is then obvious for q = 0. 
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Now, using the properties of L, we get: 



e T (i,q) > jf V s Lr-^)ds 

> rv L mus 



' t 



V 



> (T — t)VL ' 



' t 



> ¥_{T-t)L 



ff,g( s ) ds 
V T-i 

Q 



y(T-i), 

Hence, if q ^ 0, the superlinearity of L gives lim^ T Brit, q) = +00. □ 



The singularity exhibited in Proposition HI along with the limited regularity of 
L, make the use of classical theorems impossible as far as the regularity of 9t is 
concerned. We shall prove using tools from convex analysis and results on semi- 
concave functions (with nonlinear modulus of semi-concavity) that 9t is in fact C 1 
on [0, T) x R and solves a Hamilton- Jacobi equation in the classical sense. 

As a first step, we prove that 9? is locally Lipschitz on [0,T) x R and therefore 
continuous on [0, T) x R. 

In order to do that, let us prove a lemma that provides a bound to the gradient of q~ „. 



Lemma 1. 

VQ > 0, Ve G (0,T),3C(Q,e) > 0,V(t,g) G [0,T-e}x[-Q,Q], sup \q* U (t)\ < C(Q,e) 

te[t,T] 

Proof: 

We consider the case q > 0. The case of a negative q works exactly the same. 

Using the same method as for Proposition [31 we know that q- . solves an ODE of 
the form: 

p(t) = 1 ^q(t) 



where p is a C 1 function. 

From the ODE, we know that Vt G [i,T),p(t) G \p(i),p(t) + ja 2 {T - i)q). This gives 
Vte[t,T],q* u (t) G \VH>(p(t)),mm(0,VH'(p(i) + 10 - 2 (T -t)q))}. 
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Now, if H'(p(i)+-fa 2 (T -i)q) > 0, then H'(p(i) + ^a 2 QT) > and p(i) is bounded 
from below by a constant that depends on Q. 

Otherwise, we have: 



~q = ? iA (t)dt < VH'(p(t) + 7<r 2 (T - t)q)(T - t) < VH'(p(t) + ^a 2 QT)(T - t). 

Hence H'(p(i) + ^a 2 QT) > — v( £_^ > — ^ and p(t) is bounded from below by a 
constant that depends on Q and e. 

Now, in general, because Vt G [t,T],q~ Jt) > VH'(p(i)), we have a lower bound for 
qfjt). being an upper bound, we have proved the result. □ 

t,q 

Given Q > and e G (0, T), this lemma and the results of Section 2 allow to redefine 
6 T (t, q) for (t, g) G [0, T - e] x [-Q, Q] by: 



«ec 1 ([t,ri),9(f)=? )9 (r)=o ) |5(t)|<c(Q, e )7 f \ V V s J 2 / 

Using this new definition, we can prove that 9 T is locally Lipschitz on [0, T) x R: 

Proposition 7 (Local Lipschitz property). WQ > 0, Ve G (0, T), ^ Lipschitz on 
[0,T-e] x [-Q,Q]. 

Proof: 

Let us consider (t, g) and (t, g) in [0, T — e] x [— Q, Q]. 
We have: 

Denoting = + frjf ~ ?)> we nave: 



*r(t, 9) < / ( ( * — ^ + U/s) + — i (q-q)) ) da. 



Hence: 



0t& q) - ftr(f, q) < f T V s ( L ( ^) ~ L ) ) ds 
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Because L is C 1 , it is locally Lipschitz and there exists therefore a constant K(Q, e) 
such that: 

Vs [t ' T] ' L ^ — ^ — J " L \~vr ) - KiQ ' e) vj-j- 

Similarly, there exists a constant K(Q) such that: 

V* G [t, T], (q^s) + ^Z£(qf - q)j - q^sf < K(Q)\q - q\. 
Hence, there exists a constant K'(Q,e) such that 

9 T (i,q)-9 T (t,q)<K'(Q,e)\q-q\T. 
Reversing the roles of q and q we eventually obtain that: 



\9 T (i,q)-e T (t,q)\<K'(Q,e)\q-q\T. 



Now, if t < t: 



e T (t, q)= I I v s l ( ^ j + ± 7 <*$«(*) 2 j <fc + M*~ (*))• 

Since . and 5? are uniformly bounded, there exists a constant M(Q, e) such that: 



T (t,g) < M(Q,e)(t-i)+e T (t,qlfi)) < M(Q,e)(t-i)+9 T (lq)+K\Q,e)\q^(t)-q\T 

< M(Q, e)(t - t) + 9 T (t, q) + K'(Q, e)C(Q, e) |t - £|T. 
Hence, there exists a constant K"(Q,e) such that: 

< 9 T (t,q)-9 T (t,q) < K"(Q,e)\t-t\. 
Combining the results we obtain the Lipschitz property. □ 

Now, we exhibit the Hamilton- Jacobi equation solved by 9 T : 

Proposition 8 (Hamilton- Jacobi equation). 9 T is a viscosity solution of the Hamilton- 
Jacobi equation: 

-d t 9 T (t, q) - ^a 2 q 2 + V t H(d q 9 T (t, q)) = 0, on [0, T) x R. 
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Proof: 



We start with the fact that 9t is a subsolution of the Hamilton- Jacobi equation on 
(0,T) x R. 

Let (t,q) G (0,T) x R and let tp G ^((O^T) x R) be such that 9 T - <p has a 
local maximum in (t,q). We can then find r > and 77 > such that W(t',q') G 
[f - r, i + r] x [g - 77, g + 77] : 

e T (t,q)-<p(t,q)>9 T (!r,<fi-<p(!f,<fi. 

Now, let us consider i> G R. If h G (0, r] is small enough, we have that t + h < T 
and q — vh <E [q — r], q + 77]. Then: 

ip(t + h,q — vh) — ip(t,q) > 9 T (t + h,q — vh) — 9 T (t,q) 

Hence if h — > 0, we get: 

dt<p{t,q)-vd q <p(t,q) > -V t L (^j - \l^ 2 q 2 . 

Hence, Vv G R, 

-d tl p(t, q) - )p(J 2 q 2 + vd ql p(t, q) - V t L (j^j < 0. 

This gives: 

-dt(p(t, q) - \lo 2 q 2 + sup vd q (p(t, q) - V t L ( ^ ) < 

1 f V \ 

-dt<p(t, q) - ~70"V + V t sup Tydgip^ q)- L \T7) - 

-%>(*, ?) - ^¥ + V t if g)) < 0. 
This proves that 9 T is a subsolution of the Hamilton- Jacobi equation on (0,T) x R. 

We now prove that 6t is a supersolution of the Hamilton- Jacobi equation on (0, T) x 
R. 

Let (t,q) G (0,T) x R and let tp G C\(Q,T) x R) be such that 9 T - cp has a 
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local minimum in (t, q). We can then find r > and 77 > such that V(t',g') G 
[t - r, t + t\ x [q - 77, q + 77] : 

e T {t,q)- ( p{t,q)<e T {t',q')-cp{t',q'). 

Ve > 0, V7i G (0,min(r, ^)), 3(v e s ' h ) se[tjT ] G A),det0,g) such that: 

9 T (t, q) -9 T (t+ h, q£ h h ) > f t +k (v s L f^j + ^ 2 (grY ) ds - eh, 
where qf' = q — J". s v^dw. 

Now, because of Lemma [T], we can always suppose that \v e,h \ is uniformly bounded, 
independently of e and h, by a constant M. Hence, if we have h < r, h < ^=£ ? and 
Mh < 77, then G [q — 77, g + 77] and: 

?) - <P (t + h, q£l) > (v s L (^Pj + \lv 2 {qT h ?) ds - eh. 

This gives: 
Hence: 

f 1 (-fy (s,qf h ) - \ia\qty + vfd q v (s,qf h ) - V S L ds > -eh 

[~d tV (s,qr) - \lo\qfy + V S H (d& (sX^))) ds > -eh 

\^ (-dep (s,qf h ) - l^qfy + V S H (d q p (*,<£"*))) ds > -e. 

Using, the continuity of d t <p, d q (p, and H, and the uniform bound on v e,h , we then 
obtain, sending h to 0, that: 

-dt<p {t, q) - ]p<J 2 q 2 + V t H {d q p (t, q)) > -e. 
This being true for any e > 0, we have that: 

-dt(p (t, q) - )po 2 q 2 + V t H (d q <p (t, q)) > 

This proves that 9t is a supersolution of the Hamilton- Jacobi equation on (0, T) x R. 
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Consequently, 9t is a viscosity solution of the Hamilton- Jacobi equation on (0, T) xl. 
We need to prove now that the result holds on [0, T) x R. 

Let gel and let y? G C 1 ([0, T) x R) be such that 9 T — (p has a strict local maximum 
in (0, g). Then, for e > sufficiently small, (t, g) e(0,T)xI^ 9 T (t, q) - tp(t, q) - f 
has a local maximum in (t e , q e ) G (0, T) x R, and lim e _ >0 (^, Qe) = (0, g). Using what 
we already proved, we obtain: 

-dt<p(t e , Qe) + ^- \l<J 2 q\ + V t H (d q <p(t e , q e )) < 0. 

Hence: 

-dt<p(t e , q e ) - -W 2 Qe + V U H (dMte, Qe)) < 0, 

and sending e to we obtain: 

-dt<p(0, q) - ^a 2 q 2 + V H (c^(0, g)) < 0. 
Therefore, 9t is a subsolution of the Hamilton-Jacobi equation on [0,T) x R. 

Similarly, let g G R and let y? G C^QO, T) x R) be such that 9 T — (p has a strict local 
minimum in (0, q). Then, for e > sufficiently small, (t, q) G (0,T)xIh <?) - 
V?(*> ?) + f nas a local minimum in (t e ,q e ) G (0,T) x R, and rim e _>o(£e, g e ) = (0, g). 
Using what we already proved we obtain: 

-dt<p(t e , Qe)-^~ \lo 2 q\ + V t€ H (dMte, Qe)) > 0. 

Hence: 

-$¥>(t e , g e ) - -w 2 q 2 + U tE iJ (c^(* e , g e )) > 0, 
and sending e to we obtain: 

-dMO, q) - ]pa 2 q 2 + V H (d q <p(0, q)) > 0. 
Therefore, 9t is a supersolution of the Hamilton-Jacobi equation on [0,T) x R. □ 

We are now going to obtain convexity and semi-concavity properties for 9r(t, •), 
Vf G [0,T). 

Proposition 9 (Convexity of 9x(t, •)). Vt G [0, T), ^ T (t, •) a convex function. 
Proof: 

Let us consider g G R and h > 0. 

We consider q = \q* z - . + ^g? - , . and we have: 

9 T (t, q + h)- 29 T (t, q) + T (i, g - /i) 
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+ h a2 l ( q i^ 2 - 2q ^ 2 + q h-^ 2 ) ds - 



Using the convexity of L and the convexity of x — > x 2 , we obtain that: 



6 T (i, q + h)- 26 T (i, q) + 9 T (t, q - h) > 0. 
Hence, #t(£, •) is a convex function. 



□ 



Proposition 10 (Local semi- concavity of 9x(t, •)). Vt G [0, T), #t(£, •) is a locally 
semi-concave function. 



Proof: 



Let us consider Q > 0, q± < q 2 G [— Q, Q] and A G (0, 1). 



Let us define q x = ?|^ 1+(1 _ A) ^-(1-A)y (g 2 -gi) and g 2 = g 
Then: 



A0 T (*, 9i) + (1 - A)0 r (i, 92) - 0r(*, Aft + (1 - A)g 2 ) 

< / K|al('^V(i-a)l^Vl^ W(1 - a) - ,s) 



+ 



1 /" T 

27^ 2 y_ (Agi( S ) 2 + (1 - A)g 2 ( S ) 2 - q* iMl+{1 _ x)(j2 (s) 2 ) ds. 



Because L is C 1 , it is locally semi-concave and there exists a continuous and increas- 
ing function w with w(0) = such that (for e = T — i): 



Vpi,p 2 G 



-£ (C(Q, e) + 2?) , £ (C(Q, e) + 2?) , VA G (0, 1), 



XL ( Pl ) + (1 - A)L (p 2 ) - L (A Pl + (1 - A)p 2 ) < A(l - A)| Pl - p 2 \cu(\ Pl - p 2 \). 



Hence: 



X9 T (t, q x ) + (1 - X)9 T (t, q 2 ) - 9 T (t, X Ql + (1 - X)q 2 ) 



< A(l - A) J V s 



Q2(s) qi(s) 



V s V s 



u> 



92 (s) gi(s) 



ds 
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+\ia 2 \(l-\)(q 2 - qi ) 2 [ (^—^) 2 ds. 



t 



xe T {t, + (i - x)e T {t, q 2 ) - e T {t, \ Ql + (i - a)<? 2 ) 

<?2 - qi\\ , 1_ _ 2 w 1 xW . ^2 



Consequently: 



< A(l - A) |g 2 - qi\ u y — —y — J + -10- A(l - A)e(g 2 - <?i 

This proves the semi-concavity of 6*t(£, •) on [-Q, Q}. This result being true for all 
Q, the local semi- concavity of 9r(t, •) is proved. □ 

Now, we prove an important result of this section stating that 9t is in fact a C 1 
function on [0,T) x R solving the PDE of Proposition [H] in the classical sense. 



Theorem 5 (Regularity of 9t and Hamilton- Jacobi equation). 9 T £ C 1 ([0,T) x 
and: 



V(t, g) e [0, T) x E, -d t 9 T (t, q) - ]pa 2 q 2 + V t H(d q 9 T (t, q)) = 0. 



Proof: 



A first step in the proof consists in recalling that a function which is both locally 
semi-concave and convex is C 1 (see for instance Theorem 3.3.7 of |9J). This gives 
Vt G [O,T),0 T (t,.) G C^R). 

9 T is C° on [0,T) x R with Vt G [0, T), 9 T (t, ■) convex and C 1 . Hence, we deduce 
from Theorem 25.7 of [22J that d x 9j- is in fact a continuous function on [0,T) x M. 

Now, let us consider q G R. We know from Proposition [7] that 9r{-,q) G ^^(O, T), 
and then it is almost everywhere differentiable on (0, T). If t is such that 0t("j<z) is 
different iable at £, then 6*t is differentiable at (t, g) and we know that —d t 9 T (t, q) — 
\lcF 2 q 2 + V t H(d q 9T(t,q)) = 0. But this gives that almost everywhere on [0,T), 
dt9r(-,q) is equal to a continuous function. Therefore 9T(-,q) is in fact C 1 ([0,T)) 
with d t 0T(t,q) = —\lo- 2 q 2 + V t H(d q 9 T (t, q)). We conclude that <9i#r is a contin- 
uous function. Therefore 9 T G C 1 ([0,T) x R) and 9 T solves the above PDE on 
[0,T)xR. □ 

We proved that the function 9t is a smooth function on [0, T) x R and we exhibited 
the Hamilton- Jacobi equation solved by 9t- We now go further in the study of 9t 
and we start with monotonicity properties of 9t- 



Proposition 11. VT > 0,Vg G R, ^t("j?) is an increasing function on [0,T). 
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Proof: 

Consider t, i G [0, T) with t <t. Let us define q : s G [t,T] h-> 9? (s - t + *)l s < T+t . 
We have: 

< f + " fa m) + * 



* r(^(¥) + ^ (s)2 )^ 

Hence, T (-,q) is an increasing function on [0, T). □ 

Corollary 4. Consider t G M+ and g £ 1. T/ien T G (t, +00) >->■ 0T(t,q) is a 
decreasing function. 

Proof: 

Consider T, V G (t, +00) with T < T". We have that 9 T (t, q) = 6 T ,(T' -T + t,q)> 
e T '(t,q). □ 

These monotonicity results with respect to time will be useful to study the asymp- 
totic behavior of the value function. We now turn to monotonicity results with 
respect to q. 

Proposition 12. VT > 0, W G [0,T),q 6l + 4 T {i,q) is an increasing function. 
Proof: 

Consider q, q G M + with q < q. Then: 
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Mta < /(^(l¥) + ^ 2 (f^)> s 

< 9 T (t,q). 



Hence, 9 T (i, •) is an increasing function on R + . □ 

A natural corollary of these monotonicity results concerns the dependence of the 
liquidity premium on T and q: 

Corollary 5. For q > 0, the liquidity premium qS — P(T, q, S) = ^q 2 + ^q + 9 T (0,q) 
is an increasing function of q and a decreasing function ofT. 

Regarding the price of a block trade, in addition to the above natural monotonicity 
properties, two questions arise: 

• Since an agent buying a portfolio with q shares may not be time-constrained 
to liquidate it, we have to study the limit of P(T, q, S) as T — >■ +00 - this 
limit being then an upper bound for any other price with T finite. 

• Since 9t has a singularity at time t — T, it cannot be approximated using 
a numerical scheme on the Hamilton- Jacobi PDE solved by 9 T - At least not 
directly. The question of the numerical approximation of 9 T and P(T, q, S) 
requires to approximate 9t by another function without singularity. 

4.2 Block trade pricing in the time-unconstrained case 

So far, we only defined the price of a block trade when the agent buying the port- 
folio was constrained to liquidate it on the market within a certain time T. It is 
also interesting to focus on the time-unconstrained case. To that purpose, one needs 
to study the asymptotic behavior of P(T, q, S) as T — > +00, or equivalently the 
asymptotic behavior of T h-> 9t- 



Theorem 6 (Asymptotic behavior of $t). Let us consider the case of a constant 
market volume curve V t = V . 

Vg > 0,Vf > 0, T hmJ T (t,q) = O^q) = f H~* (^') dx, 
where H~ l is the inverse of H : M + — > M + . 
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Proof: 



For (t, q) G R+ x R + , T > £ i-> #r(£, 9) is a decreasing function bounded from below 
by 0. Hence linij-^+oo #r(i, g) exists. 

Since for < t' < t, Oxit 1 ,q) = 9 T -t>+t(t,q), limr-^+oo 9r(t, q) is in fact independent 
of t and we define 9 QO (q) = lim^+oo 9r{t, q). 

Using Proposition [HI ^ is a convex function and it is therefore continuous. 

Hence, for a fixed t > 0, using Dini's theorem, the convergence of (#t(£, -))t towards 
6*00 is locally uniform on R + . Now, because of Proposition [TT], we have that {6t)t 
converges towards 9^ locally uniformly on R + x R + . 

Now, we are going to use the Hamilton- Jacobi equation of Proposition [H] to prove 
that 9oo is a viscosity solution of: 

-^o- 2 q 2 + VH(9' oo (q)) = 0, q > 0, MO) = (*). 

We consider q > and (p G C 1 (R^) such that 9^ — (p has a local strict maximum in 
q. Then, we consider a sequence of triplets (T n , t n , q n ) G R+ x R+ with t n < T n such 
that lim n T n = +oo, (t n ,q n ) n converges with lim„g n = q, and 9t„ — <p has a local 
maximum in (t n ,q n ). We therefore obtain that: 

~wV n + VH(p'(q n )) < 0. 
Hence, considering the limit n — > +oo, we have: 

_I 7(T y + < o. 

Conversely, if we consider q > and G C 1 (R^) such that 9^ — tp has a local 
minimum in q. Then, we consider a sequence of triplets (T n , t n , q n ) G R+ x R^ with 
t n < T n such that lim„T„ = +oo, (t n , q n ) n converges with lim n q n = q, and 9t„ — <p 
has a local minimum in (t n , q n ). We therefore obtain that: 

+ VH(<p'(q n )) > 0. 
Hence, considering the limit n — » +oo, we have: 

-\l° 2 q 2 + VH{<p'{q)) > 0. 
This, and the fact that 9 O0 {Q) = 0, proves that 6*00 is indeed a viscosity solution of 
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the equation (*). 



Now, 6*00 is convex and hence we have almost everywhere that — ^a 2 q 2 +VH{9' 00 {q)) = 
0. Since 6^, is an increasing function (as a limit of such functions), we can write 
that almost everywhere: 



where if 1 is the inverse of H : R+ — > R + , as defined in the above statement. 



Since 9^ is convex we can integrate this equation to obtain: 



x 2 )dx = H 1 ( %—x 2 ) dx 



2V J J \2V 
and the result is proved. □ 



A straightforward consequence of this result is the following Theorem that gives the 
price of a block trade in closed-form in the time-unconstrained case: 

Theorem 7 (Closed-form expression). Let us consider the case of a constant market 
volume curve Vt — V. 

The price of a block trade for a portfolio with q > shares is: 

P(q, S) := lim P(T, q, S) = qS - \q 2 - # - f H~ l (^x 2 ) dx, 
where if -1 is the inverse of H : M. + — > M + . 

As an illustration, we can compute the price of a block trade when the execution 
cost function L is L(p) = ^|p| 1+< ^: 

Proposition 13 (Block trade pricing for L(p) = r]\p\ 1+ ^). Let us suppose that 
Vt — V and that L(p) = r]\p 
Then: 



1+0 



P(q, S) = qS--q - tpq — - — 7777 q 1+ * 



2 rhih 1 + 30 V2^ 



Proof: 



When Lip) = r/\p\ 1+ ^, the function H is given by Hip) = -r|p| 1+0 ■ Hence, 

(1+0) <t> r)<t> 



for x > 0, H^ 1 (x) = rj 1 +<i>^f-x 1 +'t> . 
Using now the result of Theorem [7J we obtain the formula. □ 
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4.3 Approximation of 9t 

In order to approximate the function $t, we are going to consider a control problem 
slightly different from the initial one in which the end point was constrained to be 
equal to 0. Now, we consider trajectories for the inventory that do not necessarily 
correspond to liquidation strategies. However, we penalize large inventories at the 
final time T. 

Given K > 0, we define for (i, q) G [0, T] x R: 

9 K , T (t, q) = inf / (v s L + ^aVis)) ds + Kq(T) 2 . 

where W^' 1+t (t 1 ,t 2 ) is the set of functions q G W 1 ' 1+e (t 1 ,t 2 ) satisfying q(ti) = a. 

Using the same ideas as above (or sometimes even more classical methods since there 
is no singularity at time t — T), we can obtain the following theorem: 

Theorem 8 (The relaxed problem). Let us fix K > 0. V(f, q) E [0,T] x R, the 
infimum that defines Oxrii, q) is attained by a unique function q* f . G W~' 1+e (t, T) 
that is in fact in C ([0, T\). 

Coming to the value function, 9k,t is a C 1 ([0, T] x R) function, solution of the PDE: 

V(f , q) G [0, T] x R, -d t 9 K , T (t, q) - ]pa 2 q 2 + V t H(d g 9 K , T (t, q)) = 0, 
with terminal condition 9x,T{T,q) = Kq 2 and "boundary" condition 9x,T(t, 0) = 0. 

To obtain a differential characterization of 9k,t, the Hamilton- Jacobi equation can 
be written in another form that guarantees uniqueness: 

Proposition 14 (Another Hamilton- Jacobi equation). Consider Q > 0. 9k,t is the 
unique continuous viscosity solution of: 

-dt9 K , T (t, q) - ]pa 2 q 2 + V t H {(d q 9 K , T (t, q)) + ) = 0, on [0, T) x (0, Q}. 

9k,t(W) = 0,Vt G [0,T], 9 K , T (T,q) = Kq 2 ,Vq G [0,Q]. 

Proof- 
Let (t,q) G [0,T) x (0,Q] and up G ^([O,^ x (0,Q]) be such that k ,t ~ <P has a 
local maximum on [0, T) x (0,(5] in (t,q). Then: 

d t 9 K ,T( t ' 0) ~ d Mt, q) < 0, d q 9 KtT (t, q) - d q ip(t, q) > 0. 
Now, since p H{p + ) is increasing, we have: 
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-d tV (t, q) - - 7 a 2 g 2 + V t H((dMt, q))+) 

< -d t e K At,q) - \i°\ 2 + v t H((d q e K , T (t, q )) + ) = o. 

This proves the subsolution part. 

Let now (t, q) G [0, T) x (0, Q\ and (p G ^([O, T) x (0, Q\) be such that k ,t - V has 
a local minimum on [0,T) x (0, Q] in (t,g). Then: 

dtOxAt^ l) ~ d Mt, V) > 0, d q e KtT (t, g) - d q ip(t, q) < 0. 
Now, since p >->■ H{p + ) is increasing, we have: 

-^(t, g) - + V t if ?))+) 

> -dte K , T (t,q) - ]pa 2 q 2 + V t H((d q 9 K , T (t, q)) + ) = 0. 
This proves the supersolution part. 

Since 6 KtT (t,0) = 0,Vt G [0,T] and 6 KtT (T,q) = Kq 2 ^q G [0,Q], # X>T is indeed a 
viscosity solution of the equation, as announced. 

Now, to prove that there is a unique continuous viscosity solution of the equation 
we use the classical doubling of variables method. 

Let us consider two continuous viscosity solutions of the equation: 9 1 and 9 2 . We 
consider a > and for e > 0, v > 0, 

U e >" : [0,T] 2 x [0,Q] 2 -> R 

(ti, t 2 , ft, ft) H- ^(ti, gi) - e 2 (t 2 , q 2 ) - - &z32)l - afa + f 2 ). 

Ve > 0,Vi/ > 0, we consider (t*'", t^, g^, 92 ") e [°; T ] 2 x [°><3] 2 such that the maxi- 
mum of U e ' u is attained in (t^, f 2 ' v , q{ v , q 2 v )- 

We have U^{?{\ t{\ q € {\ q € { v ) < U'> v (t\ v , t%", q{\ q e /). Hence: 

2^(?i - 92 ) + 2^(*i - *2 ) < # (*i > 9i ) - (*2 , 92 ) - <*(*2 - *1 )■ 

This gives that —t 2 = 0(y/e), uniformly in u, and q{ v — q 2 v = uniformly 
in e. 

Let us define B = ({T} x [0,Q]) U ([0,T] x {0}). We have to distinguish two cases: 
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• Case (i): There exists a sequence (e„,z/ n ) n converging towards (0,0) such that 
Vn G N, {t^' Vn ,q\ n ' Vn ) G B or (f 2 n ' Vn , q e 2 n '" n ) G B. 

If (t?' Vn ,ql n ' Vn ) G S then: 

U \ L 1 1 L 2 1 HI lH.2 ) 

S V [t 1 ,q 1 )—v\t 2 ,q 2 ) = V {t 1 ,q 1 ) - V [t 2 ,q 2 ). 
Similarly, if (f 2 ' Un , q 2 n ' Un ) G B then: 

TJ€n,V n (+tn,Vn Ai"n £n,V„ e n ,l/ n \ 
u \ L 1 i l 2 iH\ )H2 I 

^ y l*i >9i J _ u ih ><?2 ) — u (*i >9i > - U \ t 2 ><?2 J- 
We conclude that lim inf „^ +00 C/ £ »-""(^" , ' / ", i*"-"", gj"'"", g^'"") < 0. 

• Case (ii): For e and v below a certain threshold, we always have {t e { v , q\ v ) G" B 

and (tr , ?r) * 

In that case, we use the fact the 1 is a subsolution to obtain: 

a + - e (tr - *n - \^\r 2 + v tr H {^-{ Q r - ) < o. 

Similarly, we use the fact the 6* 2 is a supersolution to obtain: 

-a + - e (tr - *n - \i^r 2 + ^ ( 0(?r - ?n) ) > o. 

Hence, 

2« - ^ 2 (?r 2 - ?r 2 ) + - V tr) H ((^ - 92")) ) < 0. 

Sending e to and then sending z/ to 0, we get that 2a < 0. This is not possible. 

We then conclude that 

sup 9\t,q) - 9 2 (t,q) - 2at < liminf U e ' v {?{\ ?{\ q{\ q{ v ) < 0. 

(t,q)e[0,T]x[0,Q] 

Thus 9 1 - 9 2 < 2aT and sending a to 0, we get 9 l < 2 

Reversing the role played by 9 1 and 9 2 , we eventually obtain 9 1 = 9 2 . 
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Therefore, there is a unique continuous viscosity solution and the result is proved. 

□ 

The above characterization is useful to obtain numerical schemes but the main result 
of this subsection is the convergence of 9 KtT toward 9 T . This convergence result which 
is stated in the next theorem will allow to approximate 6t numerically since 9k,t can 
be computed numerically using a monotone scheme for the above Hamilton- Jacobi 
equation. 

Theorem 9 (Convergence of 9k,t)- 

V(t, q) G [0, T) x R, lim 9 K , T (t, q) = 9 T (t, q). 

Proof: 

Let us consider (i,q) G [0,T) xl By definition, 9 K ,T(i,q) < 9 T (i,q). 

In particular, Kq* f .(T) 2 < 9 T (i, q), and this leads to limx^+oo q* w t ~{T) = 0. 

Now, 

jT ^ ' g ^! )|1+e - Bv)j ds < V S L (^P^J ds < 9 K , T (t, q) < 9 T (t, q). 

Hence sup KeN II^-Li+^t) < +°o- Now, since \q* K{ J < q, we know that {q* K ^K&i 
is a bounded sequence of W~' 1+e (i, T). Hence, there exists q G W-' 1+e (t, T) such that, 
up to a subsequence, q* K $~ — > q in C°([0, T]), g(0) = 0, and <Z^ t -- ^ g in L 1+e (i,T). 

Now, using the convexity of L as in the proof of Theorem 1, we have that: 




< liminf 9 KiT (i,q) 
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Now, since 0K,T(i,q) is increasing in K, we indeed have that linift:_ H _ 00 9K,r(t : q) = 
9 T (t,q). □ 



5 Limitation to deterministic strategies: another 
route 

In Section 3, we proved that the deterministic strategy v* is in fact optimal in a 
wider class of stochastic strategies. The purpose of this section is to provide another 
proof of this result using a relaxed problem in which liquidation is not enforced. 

To start with, we introduce the set of admissible strategies of our relaxed problem: 

•4(*> q) — S (v s )se[t,T\ eV(t,T),(s,uj)^ J v T (u)dre L°°([t,T] x ft) j . 
Similarly, let us define the set of admissible liquidation strategies: 

<Ao(t,q) = j(v s ) s6 [ t ,T] e A(t,q), jf v s ds = gj . 

Proposition 15 (Verification theorem for the relaxed problem). Let us consider 
K>0. 

Let us define for (t, x, q, S) G[0,T]xlx R + x R; 

u K (t, x, q,S) = - exp ^-7 (^x + qS - ^q 2 -^q- 9 KtT (t, qfj ^ . 

Then 

u K (t,x,q, S) = 



exp (x!f' v + q^Sff'" - ^q% q > v2 - ^ q ' v - Kq^^j ^ 
where Vr e [t, T] : 



sup E 

veA(t,q) 



ql' q ' v = Q-J v * ds > 
S t T > s >" = S + J adW s -kv s ds, 



t'.s 



X^ v = x+ I [v s S^ v \ /. ( - ) - )■>> 

Proof: 
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Let us consider v G A(t, q) such that V S L ^y-^ ds < +00. Using Ito's Lemma, 

and removing the superscripts to ease readability, we have: 

u K (T,X T ,q T ,S T ) 

= u K (t,x,q,S) + d t u K (s,X s ,q s ,S s )ds + ^a 2 d 2 ss u K (s, X s , q s , S s )ds 

+ (v s S s - V S L ^y-^j - ijv s ^j d x u K (s, X s , q s , S s )ds - jf v s d q u K (s, X s , q s , S s )ds 

- / kv s d s u K (s,X s ,q s ,S s )ds + / d s u K (s, X s , q s , S s )adW s 



Hence, using the expression for %, we get: 



u K (T, X T , q T , S T ) = u K (t, x, q, S)+ Ju K (s, X s , q s , S s ) ^d t 9 KtT (s, q s ) + ^a 2 q 2 ^j ds 

- Ju K (s, X s , q s , S s ) (v s S s - V S L (j^j - ip\v 3 \ - v s (S s - kq s - ip - d q 9 K>T (s, q s )) 

+ £ 'jkv s q s u K (s,X s ,q s ,S s )ds - J ^aq s u K {s,X s ,q s ,S s )dW s . 
Simplifying this expression, and using the fact that — ^(M — v) < 0, we get: 

u K (T, X T , q T , S T ) < u K (t, x, q,S)- J -faq s u K (s, X s , q s , S s )dW s , 

- lu K (s, X s , q s , S 3 ) ^-d t 9 KtT (s, q s ) - ]pcr 2 q 2 s + d q 9 K<T (s, q s )v s - V S L {^j ^ ds - 
Using the PDE solved by 9 k ,t, we obtain that 

-d t 9 K>T (s, q s ) - ]pa 2 q 2 s + d q 6 K , T (s, q s )v s - V S L (^j < 0. 

Therefore: 

u K (T, X T , q T , S T ) < u K (t, x, q,S) - J ^aq s u K {s, X s , q s , S s )dW s . 

Now, the definition of A(t, q) guarantees that (q s ) s£ \ tt T] is bounded and hence that the 
local martingale 'yaq s UK(s,X s ,q s ,S s )dW s is in fact a martingale. Consequently: 

E [u K (T, X T , q T , S T )} < u K (t, x, q, S), 

i.e.: 
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E 



exp ( - - ( X T + q T S T - -q% - 4>qr - Kq? 



< u K (t,x,q, S). 



Now 
u K (t,x,q, S) 



exp ( -7 ( x + qS - ^q 2 - if>q - 8 K)T (t, q) 
exp ( — ~ ( x + qS — —q 2 — ipq 



V,L 



gWf 



') + ^«W*) a ) ^ - Kq* Ktq {T)' 



E 



exp ( -7 ( X£ + q^Sr ~ ifiT _ V^t _ Kq? 2 



where (X*, q*, S*) corresponds to the admissible control v* K = —q* K tA - 



Combining the results, we obtain: 



UK(t,x,q,S)= sup E 

vEA(t,q) 



- exp ( -7 ( X T + q T S T - -q T - tpq T - Kq T 



□ 



In spite of the singularity associated to the liquidation problem, we now state a 
result that gives the value function of our initial stochastic optimal control problem 
and proves that the optimal deterministic strategy of Theorem [2] is in fact optimal 
in the wider set of controls Aq. 

Theorem 10 (Verification theorem and optimality of deterministic strategies). Let 
us define for (t, x, q, S) G [0, T) x R x R + x R: 



u(t,x,q,S)= lim UK(t,x,q,S) 
K-^oo 



k 



Then: 



- exp ( -7 ( x + qS - '-q 2 -ipq- T (t, q) 



u(t } x,q,S) = sup E [— exp (— 7X^ z ' -i; )] 
veA (t,q) 

= sup E [— exp (— f yXj, x ' v ')] 

v&A 0Aet (t,q) 



and this supremum is reached by the optimal deterministic strategy v* = —q_t q - 
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Proof: 



By definition we have: 



< sup E 

veA(t,q) 



cxp 



sup E [— exp (— r fXj, x,v ^ 
veAo(t,q) 



k 



Hence: 



sup E [— exp (—jX^, x ' v )] < x, g, 5). 

veAn(t,q) 

Passing to the limit as K — > oo and using Theorem [9j we obtain that 



sup E [— exp (— 7X 7 T E ' 1 ')] < x, q, S). 

veAo{t,q) 

As above, a straightforward computation indicates then that u(t,x,q,S) is nothing 



but the expected value E 



— exp ( — 7A^ 



where = -q* G A),det- 



Hence 



■it 



(t,x,q,S) < sup E [—exp (— 7X z f : ' i; )] . 

ue^4o,det(*,9) 



Since we obviously have 



sup E [— exp (— f yX^ x ' v )] < sup E [— exp (— r )X^ x,v \^ 

v£Ao,det(t,q) veAo(t,q) 



we must have: 



u(t, x,q,S) = sup E [— exp (— 7Xy X ' 1 ')] 

veAo(t,q) 

= sup E [— exp (— 'yX^ x ' v )~\ 

veAo,det{t,Q) 



and the optimal strategy is indeed v* 



□ 



6 Numerical methods 

In this section we discuss the numerical methods that may be used to compute the 
optimal liquidation strategy and/or an approximation of the price of a block trade. 
The goal of this section is not to prove any convergence result but rather to warn 
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the reader about potential problems when using numerical methods to approximate 
the solutions of our problems. 



6.1 The flat market volume curve case 

Let us first start with the simplest case where the market volume curve is flat 
(V t = V). We proved in Theorem [3] that the optimal liquidation strategy was the 
solution of a first order ODE with maybe an unknown constant. If we use the 
notations of Theorem [31 we know that there are two cases: 

• If £ i— > r -i^ is not integrable in £ = then there exists C > such that q* is 
uniquely characterized by: 

with q*(T) = 

A consequence of this result is that q* can be approximated using a finite 
difference scheme and a shooting method. More precisely, if we consider a grid 
t = < ti < . . . < t n = nAt < . . . < t N = T, we can introduce for a given 
C > the sequence (q^) n defined by: 

q C N = 0, Mn < N , q c n _ x = q c n + VAtr' 1 (^qf + 6) 

Then, we easily see that q$ is a continuous and strictly increasing function of 
C and we can therefore easily find C numerically such that q$ = qo- 




If £ r -i 1 (g2- ) is integrable in £ = 0, then the method depends on q . If q is 
strictly larger than g im defined by: 

POO 

Aim i \f I 1 



qtfa = +oo, if / r-a— rd£ < T, 

L° t — 2 — \d£ = T, otherwise 

the same shooting method as above can be used. 

On the contrary, if q is lesser than g im , we consider the following scheme: 



g = go, Vn > 0, q n+ i 



q n - VAtr' 1 (r^qS} 



This method works very well in practice because in the case of a flat market volume 
curve we exactly know when to consider a forward scheme or a backward scheme. 
In the case of a non-constant volume curve, there is no such result and in any case 
the problem does not boil down to a first order ODE. 
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6.2 The Hamiltonian system 

In the general case of a non-constant market volume curve, we exhibited two differ- 
ential characterizations, namely the Euler-Lagrange equation and the Hamiltonian 
system of equations. We do prefer to use the Hamiltonian system to approximate 
numerically the solution. 

Let us recall that q* is solution of 



g(0) = q , q(T) = 



p(t) = 7 <7 2 g(0 
q(t) = V t H'(p(t)) 

In order to compute q* numerically, we consider the backward Cauchy problem: 



where px is a (negative) constant and we want to find px such that g(0) = go- 
Numerically this leads to the following scheme: 



\/n < N,i£-i = Pn~ 7^At^ a \ a _ n 

VnKN^qt, = q* ~ V n AtH>(p*) ~ Qn ~ 

One can easily show that for g > there exists a unique A < such that q$ = g 
(the function A < i— > q$ is continuous and strictly decreasing). However, finding 
such a A is sometimes complicated. The reason for that can be well understood in 
the case of a flat market volume curve. We know indeed in that case that some 
optimal trajectories reach before time T (this is the reason why we used a forward 
equation when q was small in the flat market volume curve case) and therefore that 
p(T) = 0. Hence A must be really small to obtain q$ = go- 
When the appropriate A seems to be unreachable numerically, a solution can be to 
consider the following forward scheme: 

Vn<N,p» +1 = ^ + 7 a 2 Atg£ 

Vn<iV,g£ +1 = (<£ + V n AtH'(j£)) + Po » % % 
and then to find numerically a minimum to the function 

iV ~ 1 / n _ n \ 1 
»<0* T N (q») = ^ V n L (\|f i ) + 

X N being a discrete counterpart of X. 
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6.3 The discrete problem 

An approach that is similar to the above one consists in considering the function 

N-l / _ \ I 

(q n )o<n<N ^ Xjv(g) = V ™ L y n vAt +1 ) + ^f 1 ^^' 
and to find a minimum to this function while imposing q = q and = 0. 

This leads to the following equation (that is obviously a discretization of the Euler- 
Lagrange equation): 

«<„<*-!, i '(^±l)- i '(£^) +V A t9 „ = 0, 5 „ = ,„, 9W = 0. 

Then, the backward shooting approach consists in introducing a sequence {qn)o<n<N 
with q v N = 0, q v N _i = v > and: 

£-i = € + Vn-^tH' (l' (^f 1 ) + . 

q% is a continuous and strictly increasing function of v and there exists a unique 
v > such that = g - However, for the same reason as above, the appropriate v 
can be really small and unreachable numerically. In such a case, we recommend to 
use the method presented above with the Hamiltonian system. 

6.4 Hamilton- Jacob i equation 

The preceding methods can be used to approximate the optimal liquidation strat- 
egy numerically. Then, if we are interested in the price of a block trade, we can 
approximate it using the function X N . It is the discrete counterpart of the two-step 
approach of Section 4. Instead, the direct approach can be used by computing nu- 
merically 6k,t for K large. 

To that purpose, we refer to the classical literature on monotone finite difference 
schemes for first order Hamilton- Jacobi equations. In our case, a scheme to approx- 
imate 9k,t on the grid {(t„,x m ),0 < n < N, < m < M} where t = < t\ < 
. . . < t n — nAt < . . . < = T and xq — < x\ < . . . < x m — mAt < . . . < xm = 
MAx = q can be: 



nn—l,m 
KST 



= n ^+--fa 2 Atx 



■V n AtH 



KST 



■\n,m- 
'K,T 



Ax 



l<n<N,l<m<M 
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0jgP = Kx 2 m , < m < M Ok°t = Q, < n < N, 

( K ' T \x T ) ) 1 - 71 - N > 1 - 

m< M. + 



Conclusion 

This article was motivated by the lack of a general framework for liquidation prob- 
lems and by the need to give a price to block trades. We developed such a gen- 
eral framework and proved existence and regularity results for optimal liquidation 
strategies. We provided two different proofs that considering stochastic liquidation 
strategies does not improve the result obtained with the best deterministic liqui- 
dation strategy. Regarding block trade pricing, we provide completely new results 
that now allow to quantify liquidity premia. Specifically, we provide a closed-form 
formula for the price a block trade when there is no time constraint to liquidate, 
and a differential characterization in the time-constrained case. 
We also complemented our results with important remarks concerning the numerical 
approximations of optimal liquidation strategies and block trade prices. 

The results we obtained can be generalized straightforwardly when one introduces 
a constant drift in the price process. Generalizing the results to the case of a multi- 
asset portfolio is currently a work in progress. 
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